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ABSTRACT 

We present a study of the stellar populations in two faint M3 1 dwarf satellites, Andromeda 
XI and Andromeda XIII. Using archival images from the Wide Field Planetary Camera 2 
(WFPC2) onboard the Hubble Space Telescope (HST), we characterize the horizontal branch 
(HB) morphologies and the RR Lyrae (RRL) populations of these two faint dwarf satellites. 
Our new template light curve fitting routine (RRFIT) has been used to detect and charac- 
terize RRL populations in both galaxies. The mean periods of RRab (RRO) stars in And XI 
and And XIII are < P„t >=0.621 + 0.026 (errorl) + 0.022 (error2), and < P„b >=0.648 + 
0.026 (errorl) ± 0.022 (error2) respectively, where "errorl" represents the standard error of 
the mean, while "error2" is based on our synthetic light curve simulations. The RRL pop- 
ulations in these galaxies show a lack of RRab stars with high amplitudes {Amp{V) > l.Q 
mag) and relatively short periods {Pab ~ 0.5 days), yet their period - V band amplitude (P- 
Amp( V)) relations track the relation defined by the M3 1 field halo RRL populations at ~ 1 1 
kpc from the center of M3 1 . The metallicities of the RRab stars are calculated via a relation- 
ship between [Fe/H], Log Pab, and Amp(V). The resultant abundances ([F e I H^Andxi - -1-75; 
[Fe/H]Andxin - -1.74) are consistent with the values calculated from the RGB slope indicat- 
ing that our measurements are not significantly affected by RRL evolutionary away from the 
zero age horizontal branch. The distance to each galaxy, based on the absolute V magnitudes 
of the RRab stars, is (m - M)o,v=24.33 + 0.05 for And XI and [m - M)o,v=24.62 + 0.05 for 
And XIII. We discuss the origins of And XI and And XIII based on a comparative analysis of 
the luminosity-metallicity (L-M) relation of Local Group dwarf galaxies. 

Key words: HB morphology, RR Lyrae, luminosity-metallicity (L-M) relation. Local Group 



1 INTRODUCTION 

Since the Sculptor and Fornax dwarf spheriodal (dSph) galaxies 
were first discovered by Shapley (1938), this class of dSph galaxies 
has drawn a significant amount of attention recently; this is because 
of their potential role in the process of galaxy formation and evolu- 
tion as the basic building blocks of giant galaxies in the cold dark 
matter (CDM) paradigm with a cosmological constant (A, White & 
Rees 1978; Hernquist & Quinn 1988, 1989; White & Frenk 1991). 
The dSph galaxies are the least luminous and appear to be the most 
common type of dwarf galaxy. It is generally believed that all dSph 
systems in the Local Group harbor ancient stellar populations as 
old as the oldest Galactic Globular Clusters (GGCs) indicating that 
Local Group galaxies may share a common epoch of early star for- 
mation (Grebel & Gallagher 2004). The total masses (dark matter 
+ baryonic mass) of the Local Group dSph systems within their 
central 300 pc are generally on the order of ~ 10^ Mq (Mateo et al. 
1993; Strigari et al. 2008; Walker et al. 2009) ,which is compara- 
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ble to the masses of the most massive Galactic GCs (i.e. the most 
massive GGC, co Cen, has a total mass of ~ 10' Mq). Despite their 
narrow range of masses, these dSph galaxies exhibit a wide vari- 
ety in their star formation histories, which was described by Mateo 
(1998) in his review as "no two Local Group dwarfs have the same 
star formation history". 

With regard to the mass range of dwarf galaxies, the recent dis- 
covery of ultra-faint dwarf (uFd) satellites around the Milky Way 
(MW) and M31 has dramatically extended the lower mass limit 
of these systems (Willman et al. 2005a,b; Zucker et al. 2006a,b; 
Belokurov et al. 2006, 2007; Sakamoto & Hasegawa 2006; Irwin 
et al 2007; Walsh, Jerjen, & Willman 2007). These newly discov- 
ered uFd systems are generally fainter than My =-8.0, significantly 
dark matter dominated (at least M/L > 100), and are some of the 
most metal-poor stellar systems ([Fe/H] < -2) found in the Local 
Group (Simon & Geha 2007). The presence of these uFd galaxies 
is predicted by cosmological simulations that seek to model the for- 
mation of the first galaxies (Ricotti & Gnedin 2005); these suggest 
that some of these faint galaxies may be "fossils" of the first galax- 
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ies in which the bulk of stars formed before the reionization of the 
Universe at z~7-10 (Bovill & Ricotti 2009). 

And XI and And XIII were first discovered by Martin et al. 
(2006, hereafter M06) along with And XII via a MegaCam sur- 
vey with the Canada-France-Hawaii Telescope (CFHT). By sum- 
ming up the flux of bright members in the upper part of the red 
giant branch (RGB), they calculated a lower limit for the absolute 
magnitude of these galaxies and obtained values in the range, -7.3 
< My < -6.4. A follow-up investigation based on Keck/DEIMOS 
spectra and Subaru/Suprime-Cam imaging was carried out by 
Collins et al. (2010, hereafter CIO). Based on their data, CIO ob- 
tained a very low metallicity ([Fe/H] — 2) for both And XI and 
And XIII. CIO also attempted to measure the distance to these faint 
satellites using both the tip of the RGB and the horizontal branch 
(HB) magnitudes. Their results locate And XI at a distance of ~760 
kpc and And XIII at a distance in the range of 760-940 kpc. 

What is particularly interesting about And XI and And XIII 
is that these galaxies appear to fill the gap between the canonical 
dSph and the newly discovered uFd populations in the Luminosity- 
Metallicity (L-M) relation of Local Group dwarf galaxies. In this 
paper, we present a comprehensive study of the stellar populations 
in these two faint M3 1 satellites using deep archival images taken 
with the Wide Field Planetary Camera 2 onboard the Hubble Space 
Telescope (WFPC2/HST). This paper is organized as follows. Sec- 
tion 2 provides a short description of the data set and photome- 
try; Section 3 describes general trends in the color-magnitude dia- 
grams (CMDs); Section 4 illustrates our RR Lyrae (RRL) detection 
method and the pulsation properties of the RRL populations found 
in each galaxy; Section 5 describes the metallicity measurements 
using two independent methods (RGB slope and RRL periods); 
Section 6 presents the distance measurement; Section 7 illustrate a 
detailed HB morphology analysis; while Section 8 discusses the L- 
M relations of our target galaxies and their possible origins.; lastly. 
Section 9 presents a summary of our results. 
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2 OBSERVATIONS AND DATA REDUCTION 

The HST/WFPC2 images of And XI and And XIII are available in 
the HST archive (program ID : GO- 11084). The central regions of 
each galaxy were imaged 16 times in the F606W (~V) filter and 
22 (And XI) / 26 (And XIII) times in the F814W (~I) filter with an 
exposure time of 1200 s. The detailed observing log is summarized 
in Table I . 

We photometered the point sources in the WFPC2 images with 
the HSTphot package (Dolphin 2000) using the identical procedure 
adopted by Yang & Sarajedini (2010). Pre-constructed point spread 
functions (PSFs) for each WFPC2 passband in the Tiny Tim PSF 
library were used to facilitate the PSF photometry. Bad pixels, cos- 
mic rays, and hot pixels were removed by using the utility software 
included in the HSTphot package before performing photometry on 
the images. Aperture corrections, which are defined by the average 
difference between the PSF photometry and aperture photometry 
with a 0.5 arcsec radius, were calculated using the default settings 
of HSTphot. The instrumental magnitudes were transformed to the 
ground- based JohnsonCousins VI system using equations included 
in HSTphot. The final list of stars from the resultant photometry in- 
cludes only the point sources flagged as good stars (Object type I) 
by the HSTphot classification. 
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Figure 1. The VI color-magnitude diagram of And XI. RRab candidates 
detected from our template light curve fitting routine (see section 4) are 
mai'ked as red open circles. The right panel illustrates the photometric com- 
pleteness indicating that our photometry is better than ~90% complete at 
the level of the horizontal branch. 
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Figure 2. Same as Figure 1. but for And XIII. 



3 COLOR-MAGNITUDE DL4GRAMS 

Color-magnitude diagrams (CMDs) of And XI and And XIII in the 
VI passbands are presented in Figures I and 2. The right panels 
of each figure illustrate the photometric completeness derived from 
HSTphot's artificial star feature (Yang & Sarajedini 2010), indicat- 
ing that the point sources in both galaxies are well photometered to 
~ I magnitude below the HB level with better than ~90% photomet- 
ric completeness. Thus, it is reasonable to assume that photometric 
incompleteness does not significantly affect the rest of our analysis, 
which is focused on the horizontal branch stars and brighter. 

Both galaxies exhibit relatively steep and narrow red giant 
branches (RGBs) compared to the more luminous (My < -8) dwarf 
satellites in the Local Group. The HBs are well defined with several 
RR Lyrae candidates (see section 4). The CMDs of both galaxies 
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Table 1. Observing Log. 



Object 


R.A. (12000) 


Dec (J2000) 


Filters 


Exp Time 


Data Sets 


HJD Range (+2 454 000) 


And XI 


00 46 20.0 


+33 47 30.0 


F606W 
F814W 


16x 1200s 
22x 1200s 


u9x701-u9x702 
u9x703-u9x705 


352.78677 - 353.26802 
350.12355 -364.38978 


And XIII 


00 51 50.87 


+33 00 18.20 


F606W 
F814W 


16x 1200s 
26x 1200s 


u9x711-u9x712 
u9x713-u9x715 


304.23383 - 305.44764 
303.76789-314.70330 



also exhibit no significant signs of young (< IGyr) main sequence 
stars or intermediate age (1~ 10 Gyr) asymptotic giant branch popu- 
lations. The overall features of the CMDs for And XI and And XIII 
are reminiscent of metal-poor Galactic GCs, indicating that these 
galaxies seem to be purely old and metal-poor stellar systems and 
have experienced relatively simple star formation histories (SFHs). 



4 RR LYRAE VARIABLES 

4.1 Detection and Characterization 

We searched for RRL stars in And XI and And XIII by employ- 
ing a newly developed template light curve fitting method dubbed 
"RRFIT (Robust RR Lyrae light curve FITing)" written in FOR- 
TRAN. This new period searching routine works in a similar fash- 
ion to FITLC (Mancone & Sarajedini 2008), which is our previ- 
ous template light curve fitting code, but it is especially designed 
for large RRLs survey programs. Unlike other optimizing meth- 
ods which use a gradient search or a brute force fitting routine 
to find the best-fit model parameters, RRFIT implements a robust 
genetic algorithm known as "PIKAIA", which efficiently finds the 
best-fit model parameters via an artificial intelligence random mesh 
searching technique. RRFIT takes full advantage of the PIKAIA 
algorithm, and as a result, significantly enhances computing speed 
(e.g. RRFIT is at least 10 times faster than FITLC) without sac- 
rificing fitting accuracy. Another advantage of using RRFIT is the 
variety of standard templates that it can employ. In addition to the 
8 RRL light curve templates employed by FITLC, we have adopted 
an additional 17 unique RRab templates from the work of Kovacs 
& Kupi (2007) into our template light curve library to achieve a 
better description of various types of RRLs. 

In order to detect RRL populations in these faint dwarf satel- 
lites, we begin by identifying all of the variable star candidates 
within a range of V magnitude (24.5< V <26.0) by calculating 
the reduced Xyi value of each star defined by the following formula 
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Any peculiar data points that differ from the mean magnitude 
by ± 3 cr were excluded from the Xvi calculation. We compiled 
the list of potential variables by selecting stars with a x\i value 
greater than 3.0. For reference, the reduced x\i values of typical 
non-variable stars at the HB level (V(HB)~25.4) of And XI and 
And XIII is less than 3.0. This variability threshold generated a 
list of 69 variable star candidates (46 in And XI; 23 in And XIII) 
from the WFPC2 images of both galaxies. Then, we ran the RRFIT 
routine on the VI time series photometry of these variable star can- 
didates in order to find the best-fit light curve parameters, which 
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Figure 3. The best-lit light curves of RRL candidates in And XI. 



include period, amplitude, and the epoch of maximum light. Af- 
ter careful eye-examination of the resultant light curves, we found 
17 RRLs (10 RRab, 1 RRc and 6 unclassified) in And XI and 9 
RRLs (8 RRab and 1 unclassified) in And XIII. The best fitting 
light curves for these RRL candidates from our RRFIT analysis 
are presented in Figure 3 and 4. The pulsation properties of the 
RRL candidates in each galaxy are listed in Table 2. The mean V 
band magnitudes of the RRab candidates from the best fitting light 
curves are < V{RR) > = 25.31 ± 0.02 and 25.49 ± 0.02 for And 
XI and And XIII, respectively. The given uncertainties represent 
the standard error of the mean. We consider < V(RR) > values as 
the best estimate of the < V(HB) > value for each galaxy. Several 
RRL candidates are noted as "unclassified" because the goodness 
of fit yielded by the RRc and RRab templates are indistinguish- 
able. These tend to be the lower amplitude variables. A detailed 
discussion of these "unclassified" RRL candidates is presented in 
the following section. 

We also attempted to search for anomalous Cepheids (ACs) 
whose pulsation periods (0.4 < P < 2 days) overlap significantly 
with the RRL period range. However, neither And XI nor And XIII 
appear to contain any ACs. Since ACs are produced by young-to- 
intermediate age (1-3 Gyr) populations, the lack of ACs in these 
galaxies supports our assessment of the overall star formation his- 
tories of these dwarf galaxies as inferred from their VI CMDs. 
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Table 2. Properties of RR Lyrae stars. 





R.A. 


Decl. 






Period 






ID 


(J2000) 


(J2000) 


<V> 
And XI 


<V>-<I> 


(days) 


A(V) 


Type 


Vn040Q 
V vjyj^yjy 


4(^ ]f\ 47 






7032 


6347 


8044 


ab 




df^ \f\ M 


4_l^ 47 1 Q 47 


AJ. 1 OOU 


5168 


4655 


4905 




V00424 


D 4^1 1 S 7^ 


4-^^ 47 0^ fift 




7256 


5610 


7979 


ab 


vnn-d7s 

V UVT-ii. J 


n 4(^ 1 S 04 
U to 1 j.yt 


1 /17 I ^; 7S 


25 295 1 


5369 


5313 


n QHVQ 
U.yV 1 V 


ab 


V00447 


46 17.30 


+33 47 09.12 


25.2636 


0.5615 


0.4290 




c? 


V02120 


46 14.97 


+33 47 06.55 


25.1377 


0.7162 


0.7230 


0.6420 


ab 


V02157 


46 15.81 


+33 46 50.74 


25.2600 


0.6107 


0.8000 


0.5667 


ab 


V02172 


46 15.16 


+33 47 12.95 


25.2149 


0.6010 


0.3852 


0.3750 


c? 


V02193 


46 16.75 


+33 46 55.58 


25.3169 


0.6885 


0.5539 


0.7725 


ab 


V02195 


46 15.63 


+33 46 57.88 


25.2808 


0.6357 


0.5080 


0.4515 


alias? 


V02200 


46 16.72 


+33 46 54.48 


25.3353 


0.6331 


0.5764 


0.7936 


ab 


V02204 


46 16.48 


+33 47 21.37 


25.3839 


0.7341 


0.6327 


0.5899 


ab 


V02219 


46 16.87 


+33 47 09.78 


25.3640 


0.5801 


0.4275 


0.6388 


c? 


V02221 


46 16.86 


+33 47 08.53 


25.3588 


0.6054 


0.3793 


0.5809 


alias? 


V02224 


46 17.24 


+33 47 13.06 


25.3926 


0.6089 


0.5757 


0.9582 


ab 


V02237 


46 17.56 


+33 46 54.70 


25.4184 


0.7088 


0.6257 


0.7728 


ab 


V04338 


46 17.04 


+33 46 49.26 


25.3653 


0.6123 


0.4070 


0.5766 


c? 


And XIII 


V00026 


51 46.00 


+33 04.38 


25.5649 


0.7299 


0.5949 


0.5677 


ab 


V00620 


5146.18 


+33 11.12 


25.4153 


0.7894 


0.6759 


0.6772 


ab 


V00629 


51 44.58 


+33 08.32 


25.4071 


0.7464 


0.6134 


0.4953 


ab 


V00636 


51 46.99 


+33 08.71 


25.4605 


0.7675 


0.8116 


0.2972 


ab 


V00651 


51 46.81 


+33 07.11 


25.5733 


0.7460 


0.5980 


0.8240 


ab 


V00652 


51 46.55 


+33 16.74 


25.5419 


0.6450 


0.6046 


0.8455 


ab 


V02129 


51 46.35 


+33 13.60 


25.4872 


0.7252 


0.6725 


0.7452 


ab 


V03451 


51 46.88 


+33 09.63 


25.4999 


0.7173 


0.6158 


0.5103 


ab 


V02147 


51 46.48 


+33 12.96 


25.3780 


0.5491 


0.4287 


0.4970 


c? 



ages, ~ 1000 synthetic light curves of RRab and RRc stars were 
generated. Periods and amphtudes were randomly assigned to each 
artificial RRL within appropriate ranges for each type (RRab : 
0.5 < P < 1.2 , 0.2 < Amp < 1.5; RRc : 0.2 < P < 0.5, 
0.2 < Amp < 0.5). Photometric errors for each data point were 
assigned using gaussian deviates centered at 0.06 mag, which is a 
typical photometric error at the level of the HB for both galaxies in 
the WFPC2 photometry. Then, we ran the RRFIT routine on these 
artificial RRL stars. 

The plots shown in the top panels of Figure 5 and Figure 6 
illustrate the difference between the input and output periods as 
a function of input period. For And XI (see Figure 5), RRc stars 
appear to experience more significant aliasing compared to RRab 
stars. Quantitatively, ~ 84 % of the input periods for RRab stars 
were accurately recovered within ± 0.05 day, while ~ 61 % of the 
RRc periods were recovered within the same ± 0.05 day period 
range. In addition, the output periods of the RRc stars tend to be 
longer than the input periods. The negative tail shown in the AP 
distribution (the left panel in the middle of Fig 5) for the artificial 
RRc stars of And XI is mainly caused by this bias in the RRc pe- 
riods. This result reveals that there can be considerable uncertainty 
in the properties of individual RRc stars in And XI. 

To investigate the degree of misidentifications in the RRc 
regime, we performed a comparative analysis between the best-fit 
light curves and the first-overtone mode light curves for those un- 
classified RRL candidates found in the previous section. Observing 
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Figure 4. Same as Figure 3, but for And XIII. 
4.2 Synthetic Light Curve Simulations 

In order to quantitatively gauge the influence of false periods 
(aliases) on our analysis, we have performed the following set of 
synthetic light curve simulations (Sarajedini et al 2009 (S09); Yang 
& Sarajedini 2010; Yang et al 2010, hereafter YIO). Applying the 
same observing windows as the And XI and And XIII WFPC2 im- 
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Figure 5. The results of our synthetic hght curve simulations for the And XI 
data. Top panels illustrate the difference ( AP=P;„ - Pour) between the output 
period and the input period as a function of the input period for two different 
pulsation modes. Lower panels (middle & bottom) show the distributions of 
the deviation of the output peiiod and amplitude from the input values. Our 
simulations reveal that the resultant periods from our analysis are likely to 
have more aliasing in the RRc regime as compared to the RRab regime. 
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Figure 6. Same as Figure 5, but for And XIII. The simulations show that 
the aliasing issue is less significant for the And XIII RRL candidates as 
compared with the And XI RRLs. 
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Figure 7. Comparisons between the RRab and RRc template fits for the 
unclassified RRL candidates. Note these unclassified RRL tend to have 
smaller amplitudes than most other RRLs in our sample. 



cadence and phase coverage are the two most important factors for 
accurate period measurements (YIO), and these factors also have a 
significant role in determining which type of light curve provides 
the best fit to the observations. Photometric errors also influence 
the appearance of the measured light curves. Large photometric er- 
rors can deform the original shape of a light curve and produce a 
misidentification of the type of variable (Yang & Sarajedini 2010). 
However, based on our simulations, the measured periods and am- 
plitudes are largely insensitive to the actual classification of the 
RRL (e.g. whether it is classified as an RRab or RRc). From Figure 
7, we see that the majority of unclassified RRL stars exhibit good 
quality fits both in the fundamental and first-overtone modes. How- 
ever, their location in the Period-Amplitude (P-A) diagram (see 
Sec. 4.3) places them in the RRc regime. Hence, we conclude that 
these unclassified RRLs are likely to be RRc stars. There are two 
exceptions to this - V02195 & V02221, which show significantly 
different solutions in the resultant periods for each type of pulsation 
mode (see Figure 7). Again, relying on their location in the P-A di- 
agram, we suspect that the RRc template is more appropriate for 
these stars and that the RRab period is the result of an alias. There- 
fore, the corrected census of RRc stars in each galaxy is A^^ = 5 and 
1 for And XI and And XIH, respectively. 

In order to estimate the errors of the individual RRL periods, 
we performed the following statistical test introduced in YIO. From 
the artificial RRL lists, we randomly draw the same number of arti- 
ficial RRL stars as our observed RRLs in each galaxy [And XI : 10 
(RRab) and 5 (RRc); And XIII : 8 (RRab) and 1 (RRc)] to calculate 
an average AP value for each sample drawn. Note that since And 
XIII only has one RRc star, we have excluded these And XIII stars 
from this analysis. We construct the distribution of < AP > val- 
ues for each pulsation mode by iterating this probability sampling 
10,000 times. Figure 8 illustrates the resultant < 1S.P > distribu- 
tions of the artificial RRL stars for And XI and And XIII along 
with the best-fit Gaussian distributions. The 1-cr errors of the best- 
fit Gaussians, as given in each panel, are good approximations to 
the standard deviations of each < AP > distribution. Therefore, 
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ID 




<AP> (days) <AP> (days) 

Figure 8. The distributions of the output period deviation for synthetic 
RRab and RRc stars. The dashed lines illustrate the best-fit Gaussians for 
each distribution. < P >a,.g and tr represent the peak and the standard devi- 
ation of the Gaussian fits, respectively. 

the standard errors (cr/ VA') of the individual RRL periods in our 
analysis are cr(Pc)=±0.038 and cr(P„4)=±0.022 days for the RRL 
candidates found in And XI. For the And XIII RRab candidates, 
we find a(Pab)=±0.022 days. The same method was applied to es- 
timate the errors in the amplitudes. We obtained cr(Amp(V)c^) = ± 
0.027 and a-{Amp(VU) = ± 0.024 mag for the And XI RRL candi- 
dates. For the And XIII RRab stars, the error value is (T{Amp(V)ab) 
= ± 0.033 mag. We found no significant biases in our determination 
of the mean periods and amplitudes. Hence, we assume that period 
aliasing does not significantly affect the results in the remainder of 
our analysis, which is focused on the mean periods and amplitudes. 

4.3 Period-Amplitude (P-A) Diagrams 

The P-A relations of the RRL populations in And XI and And XIII 
are presented in Figure 9. For comparison, we generated contoured 
P-A diagrams for the RRL stars in the spheroid of M3 1 from S09 
representing a population at Rgc ~ 5 kpc) and from Brown et al. 
2004 (B04, hereafter referred to as the Brown field) at Rgc -11 
kpc. The contour plot for the Brown field was scaled up by a factor 
of 5 in order to see its trend more clearly. The RRL stars (open 
symbols) discovered in And XI and And XIII from our analysis are 
plotted on top of these contoured P-A diagrams. Two solid lines 
representing the Oosterhoff I (Oo I) and II (Oo II) sequences for the 
Galactic Globular Clusters (Clement 2000) are also plotted. Since 
the amplitudes of the M31 RRL stars in these two previous studies 
were determined in the ACS/F606W passband, we applied an 8% 
increment to the amplitudes of each M3 1 RRL in order to convert 
the F606W amplitudes into V band amplitudes (B04; S09). 

From Figure 9, we see that the P-A relations of the RRab stars 
in And XI and And XIII agree well with the general trend of the 
stars in the Brown field, while the RRc candidates in both dSph's 
are obviously shifted toward longer periods with respect to the M3 1 
RRc stars. The mean RRab periods for And XI and And XIII are 




log P (days) log P (days) 

Figure 9. The Period-Amplitude (P-A) relations of RRL candidates (open 
circles) in And XI and And XIII are plotted on top of the contoured P- 
A relations of the M31 halo RRLs (left panels, S09 : inner halo, Rg^ ~ 
5kpc ; right panels, B04 : Brown field, Rgc ~ llkpc). The contour plot for 
the Brown field was scaled up by a factor of 5 in order to see its trend 
more clearly. Solid lines represent the fiducials of Oosterhoff I and II GGC 
(Clement 2000). 



< P„t > =0.621 ± 0.026 (errorl) ± 0.022 (error2) days and 0.648 
± 0.026 (errorl) ± 0.022 (error2) days, respectively, while for the 
Brown field, this value becomes 0.594 days (B04). The "errorl" 
value represents the standard error of the mean and the "error2" 
value is the error derived from our synthetic light curve simulations. 
The difference in the mean RRab periods between the M3 1 halo and 
these two dSph's is mainly due to the lack of short period RRab 
stars with high amplitudes in And XI and And XIII. 

To better understand the characteristics of the RRL popula- 
tions of these two dSph's, we compare the P-A relations of And 
XI and And XIII with those of four other Andromeda dSph's (And 
I & III (P05); And II (P04); And VI (P02)) in Figure 10. Accord- 
ing to this comparison, the P-A relations of the RRL stars in And 
XI and And XIII show the closest resemblance to that of And III 
among these four Andromeda dSph's. The RRab stars in And III 
track the general trend of the Brown field RRab stars (B04). The 
lack of short period RRab stars with high amplitudes and the lo- 
cations of the RRc stars in the P-A relation of And III mirror the 
characteristics of the RRL populations in And XI and And XIII. It 
is interesting to see that And III is the least luminous galaxy among 
these four dSph's and its mean metallicity (< [Fe/H] >=-1.7) is 
comparable to those of And XI and And XIII (< [Fe/H] >Andxi=- 
1.74; < [Fe/H] >Andxiii=-^-T5, see section 5). It seems reasonable 
to assert that the behavior of the P-A relations of these three dSphs 
(And III, And XI, and And XIII) and their similar global proper- 
ties could be related. From Figure 10, we also see that the P-A 
relations of And I, II, III, and VI are as diverse as their complex 
star formation histories (Mateo 1998). The broad dispersion in the 
RRab periods of each Andromeda dSph also reflects the signifi- 
cant metallicity dispersions in these galaxies supporting previous 
assertions that these Andromeda dSph galaxies have likely experi- 
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Alcock et al. (2000). The relation is given by the following equation 
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Figure 10. A comparison of the P-A diagrams for the M31 dSph galaxies. 
See text for discussion. 



enced extended and complex chemical enrichment histories during 
the course of their evolution (P02, P04, P05). 



5 METALLICITY 

M06 derived the mean metallicities of And XI and And XIII by in- 
terpolating between Padova isochrones (Girardi et al. 2004) plotted 
on the CFHT/MagaCam survey data. The estimated mean metal- 
licity values are < \FejH\ >=-1.3 and -1.4 dex for And XI and 
And XIII respectively with an error of ±0.5 dex. Most recently, 
from their Keck/DEIMOS and SUBARU/SuprimeCam study of the 
faint M31 satellites, CIO estimated metallicities of individual RGB 
stars in these two dSph's via isochrone interpolation using the Dart- 
mouth isochrones (Dotter et al. 2008) as well as Calcium II triplet 
spectra. The mean metallicities obtained by applying these two dif- 
ferent methods converged at < \FejH\ >=-2.0 ± 0.2 dex for both 
And XI and And XIII. CIO explained that the relatively higher 
mean metallicity values for these dSph's in the previous work of 
M06 were probably caused by the use of the Padova models, which 
systematically produce more metal-rich values in the MegaCam fil- 
ters than other isochrone models. 

We estimated the metallicities of And XI and And XIII via the 
period-amplitude-metallcity relationship for RRab stars derived by 



[Fe/H] = -8.85 \\-0%Pab + 0.15Amp(V)] - 2.60. 

The estimated mean metallicity of the RRab stars in And XI is 
< [Fe///] >= -1.75 ± 0.12 (errorl) + 0.13 (error2). For And XIII 
this value becomes < [Fe/H] > = -1.74 + 0.12 (errorl) ± 0.14 
(error2). Errorl represents the standard error of the mean, while 
error2 is the amount of error propagated from the determination 
of individual periods and amplitudes of the RRL candidates. Thus, 
the total errors in our metallicity determinations as described by the 
quadrature sum of errorl and error2 are + 0.18 dex for both And 
XI and And XIII. To double check the validity of our estimates, we 
independently calculated the mean metallicity of these dSphs using 
the S index (Hartwick 1968; Saviane et al. 2000). This represent 
the slope of the RGB originally defined by the line connecting two 
points of the RGB in the (V,B-V) plane - the first point is at the 
level of the HB, and the other one is at 2.5 mag brighter than the 
HB. Saviane et al. (2000) redefined this relationship in the (V, V-I) 
plane by adjusting the second point to 2.0 mag brighter than the 
HB. The metallicity-S index relationship calibrated by Saviane et 
al. (2000) is given in the following form : 

[Fe/H] = -0.24 x 5 + 0.28 ± 0.12 {rms) 

Figure 11 illustrates our measurements of the S indices for 
And XI and And XIII. The (V-I) colors of RGBs at V(HB) and 
V(HB-2.0) were calculated from a least squares fit of a quadratic 
polynomial to the RGBs of each galaxy. By using the metallicity- 
S index relationship above, we obtained [Fe/H] = -1.70 and -1.63 
dex for And XI and And XIII, respectively. The metallicity esti- 
mates for both And XI and And XIII from the RGB slope show 
excellent agreement with the RRL metallicities. It is worth men- 
tioning here that we are aware of a caveat about the use of the MA- 
CHO relationship (Alcock et al. 2000). According to Cacciari et al 
(2005), this method can underestimate the metallicity especially for 
RRLs that have significantly evolved from the zero-age horizontal 
branch. This evolutionary effect generally raises the luminosity of 
RRL stars. Then, the pulsation period of these evolved RRLs be- 
comes longer, as a result, yielding a relatively lower metallicity. 
Our metallicity estimates based on these two completely indepen- 
dent methods agree very well and difi^er by only ~ 0.1 dex in the 
worst case. Therefore we conclude that the RRL evolutionary effect 
does not negatively impact our metallicity determination using the 
MACHO relationship for And XI and And XIIL 



6 REDDENING AND DISTANCE 

Considering the brightest RGB star confirmed by their 
Keck/DEIMOS spectroscopy as a proxy for the tip of the 
RGB (TRGB) in each galaxy, CIO set distance ranges of 610-770 
kpc, and 750-940 kpc for And XI and And XIII respectively. 
For this calculation, they assumed that the brightest DEIMOS 
star is found within 0.5 mag of the TRGB. CIO also provided 
distances for And XI and And XIII estimated by utilizing their HB 
magnitudes. In this case, the distances of And XI and And XIII 
become ~ 830 kpc, and 910 kpc, respectively. 

The previous studies (M06 and CIO) used a reddening value 
of E(B-V)=0.08 taken from the Galactic extinction map of Schlegel 
et al. (1998) in the calculation of the distances for And XI and And 
XIII. To double check whether the use of this reddening value is 
reasonable, we have calculated the metallicity values for each dwarf 
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Figure 11. The illustration of the S index method for the measurement of 
mean metallicity (Saviane et al. 2000). The red solid lines represent RGB 
fiducials for each galaxy derived from a least square fit of a quadratic poly- 
nomial. The points marked as large plus symbols are the V magnitude of 
the RGB at the HB level and that of the RGB at 2.0 mag brighter than the 
HB level. 

galaxy by using this reddening in the relationship between metal- 
licity and intrinsic RGB color at the level of HB (i.e. (V - /)o,g) 
derived by Saviane et al. (2000). 

The observed RGB colors at the level of HB ((V - I)rgb) for 
And XI and And XIII are 1.108 + 0.020 and 1.061± 0.013 (see 
Table 3 in the following section) respectively. Then the adopted 
reddening value (E(B - V) = 0.08) was converted into E(V - I) 
using E(V - I) = 1.38£'(5 - V) (Tamman et al. 2003) in order to 
calculate (V-I)Q,g. The metallicity-(V-/)o,g relationship presented 
by Saviane et al. (2000) is given in the following form : 

[Fe/H] = 5.25 x (V - I)o,g - 6.52. 

It should be noted that the metallicity used for this calibra- 
tion is on the Zinn & West scale (Zinn & West 1984). The adopted 
reddening yields the resultant metallicity values of [Fe/H]/i„^ixi = 
-1.28 ± 0.11 and lFe/H]A„dxiii = -1.53 + 0.07 both of which are 
significantly more metal-rich than not only our estimates via the 
RRL period-amplitude relation and the RGB slope, but also the 
metallicity values of M06 and CIO, which used the reddening value 
of E(B- V)=0.08. Especially in the case of And XI, these metallici- 
ties are clearly not consistent with the overall properties seen in the 
CMDs and RRL populations of these galaxies. This result suggests 
the need for a more robust estimate of the reddening toward And 
XI and And XIII. 

Using the metallicity values derived from the RGB slopes of 
And XI and XIII, we calculated (V - /)oj using the metallicity- 
(V-I)o,g relation from Saviane et al. (2000). These (V-I)o.g values 
were then combined with the observed (V - I)g quantities and the 
conversion between E(V - I) and E(B - V) (Tamman et al. 2003, 
see above) to give reddening values of E(B - V)A„dxi=0-i5 and 
£(S- V)A„rfx/;/=0.11. 

Moving on now to calculate the distances, we note that the 
mean V magnitude of 11 RRab candidates found in And XI is 
< V{RR) > = 25.31 ± 0.02. By applying the extinction law, Av = 
3.1E(B-V) = 0.450, we obtain a mean intrinsic V magnitude for the 
And XI RRab stars of < Vo(RR) > = 24.86 ± 0.02. Then, the abso- 
lute V magnitude of each And XI RRab star is calculated using the 



metallicity-luminosity relationship for RRL stars, Mv(RR) = 0.23 
[Fe/H] + 0.93, from Chaboyer (1999). This yields a mean absolute 
V magnitude of< Mv(RR) > = 0.53 ± 0.05. The quoted errors rep- 
resent the quadrature sum of the standard error of the mean (0.03 
mag) and the amount of error propagated from our metallicity de- 
termination (0.04 mag). Thus, we find the mean absolute distance 
modulus for And XI to be (m - M)o = 24.33 ± 0.05. This places 
And XI at a distance of ~ 734 ±18 kpc, which is consistent with the 
distance range of 610-770 kpc obtained from the TRGB method, 
but smaller than the distance of 830 kpc estimated using the HB 
magnitude with smaller reddening value in the previous work of 
CIO. 

For And XIII, the apparent mean V magnitude of 8 RRab can- 
diates is < V(RR) > = 25.49 ± 0.02. The estimated extinction, Av 
= 0.341 yields the intrinsic V magnitude for the And XIII RRab 
stars of < Vo(RR) > = 25.49 ± 0.02. Employing once again the 
Chaboyer (1999) equation for the RRab metallicity-luminosity re- 
lation, we find < Mv(RR) > = 0.53 ± 0.05. This gives a distance 
modulus of (m - M)o = 24.62 ± 0.05 (D ~ 840 ± 18 kpc) for And 
XIII, placed in the middle of the distance range of 750 - 940 kpc 
estimated in the previous study of CIO. 



7 HORIZONTAL BRANCH MORPHOLOGIES 

We make use of a newly quantified HB index, A(V-I) (Dotter et al. 
2010, hereafter DIO), in order to characterize the HB morphologies 
of And XI and And XIII. This index is defined as the difference in 
color between the HB and the RGB at the level of the HB. The 
advantage of using this metric in the measurement of HB morphol- 
ogy is twofold. First, the A(V-I) index has less degeneracy at both 
(red & blue) ends of the HB morphology range (see Fig. 2 of DIO) 
than the conventional HB index, (B-R)/(B-l-V-l-R) (Lee et al. 1994), 
which relies on counting stars on the blue (B) and red (R) portions 
of the HB as well as the RR Lyrae variables (V). Second, the cal- 
culation of A(V-I) is relatively straightforward thereby minimizing 
the uncertainties associated with the census of RR Lyrae variables 
and the contamination of the red HB by intermediate-age RC stars. 

Figure 12 illustrates the measurement of the A(V-I) indices 
for And XI and And XIII. The mean V magnitude of the HB, 
< V(HB) >, was obtained from the mean V magnitude of RRab 
stars, < V(RR) >, calculated from the best fitting light curves 
(< V{RR) >Andx, = 25.31; < V{RR) >Andxm = 25.49, see sec- 
tion 4.1). Stars inside the blue-boxed regions were considered to 
be HB stars, while RGB stars at the level of the HB were selected 
within the red rectangles in the VI CMDs of each galaxy. The dif- 
ference between the median colors of HB and RGB stars was cal- 
culated to yield A(V-I) indices for each galaxy. The uncertainties in 
the A(V-I) indices were calculated by the quadratic sum of the l-cr 
errors in the median colors of the HB and RGB. We used a boot- 
strapping method for estimating errors in the median colors. Ten 
thousand re-samples of the HB and RGB stars were constructed by 
using random sampling with replacements from the observed HB 
and RGB data. The standard deviation of the bootstrapped medians 
represents the l-cr errors in the median colors of the HB and RGB. 
The HB parameters for each galaxy are listed in Table 3. 

The [M/H] vs A(V-I) diagram for 65 GGCs (DIO) is shown in 
Figure 13. Solid circles and open squares represent inner {Rac <8 
kpc) and outer {Roc >8 kpc) halo GCs, respectively. And XI and 
And XIII are marked as blue open stars along with 10 other Lo- 
cal Group dSph galaxies (Harbeck et al. 2001; Sarajedini et al. 
2002; Komiyama et al. 2007). The solid line is an empirical fit 
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Figure 12. These panels show the measurement of the A(V-I) index. The 
blue and red boxed regions describe the selection of HB and RGB stars, 
respectively. The horizontal line presents the level of the horizontal branch 
obtained from the mean V magnitude of RRab stars. The vertical Unes indi- 
cate the median colors of the HB and RGBs. It should be noted that RRab 
stars were not included in the calculation of the median color of the HB (see 
DIO for more detail). The dotted lines are RGB fiducials for each galaxy ob- 
tained by fitting a quadratic polynomial. 



to the inner halo GCs from DIO. We adopted the [Fe/H] values 
of And XI and And XIII from section 5 of the present paper in 
which the [Fe/H] - logPab - Amp(V) relationship from Alcock et 
al (2000) is used for the metallicity derivation. Then, these [Fe/H] 
values are translated into [M/H] by applying the following relation, 
[M/H] = [Fe/H] + logio(0.638 x lO'""'^! + 0.362) from Salaris et 
al. (1993) by assuming [Q'/Fe]=0.2 (Tolstoy et al. 2003) for both 
galaxies. For the other 8 dSph galaxies, the metal abundance val- 
ues are taken from the compilation of Grebel, Gallagher & Har- 
beck (2003); the corresponding [M/H] values are calculated with 
the same procedure used for And XI and XIII. Since the HB in- 
dices of the additional 8 dSph galaxies are in the canonical form 
of (B-R)/(B-hV-i-R), these values need to be converted into the 
new A(V-I) index. Figure 14 illustrates the correlation between (B- 
R)/(B-l-V-l-R) and A(V-I) derived by a least squares fit of a cubic 
polynomial to the data set obtained from Figure 2 of DIO. We 
have used this relation to calculate the A(V-I) values for the 8 dSph 
galaxies. To double check the validity of this relationship, we mea- 
sured (B-R)/(B-)- V-l-R) indices for And XI and And XIII (see Figure 
15) and plotted them in the (B-R)/(B-hV-i-R) vs A(V-I) relation with 
an open star symbol. For the (B-R)/(B+V+R) index calculation, the 
boundaries of the instability strip in the VI CMDs were determined 
by the blue and red edges of the RRL (V-I) colors. From Figure 
14 we see that the locations of And XI and And XIII agree well 
with the derived relationship (dashed line). We summarize the HB 
indices and other physical parameters of our sample dSph galaxies 
in Table 4. 

From Figure 13, we see that the Local Group dSph galaxies 
and the GGCs show different behavior in the metallicity-HB mor- 
phology diagram; that is to say the dSph galaxies exhibit HBs that 
are generally redder than those of the GGCs at the same metallic- 
ity. This suggests that the stellar populations in these Local Group 
dwarf galaxies are subject to the second parameter effect. If we ac- 
cept the conventional age hypothesis (Gratton et al. 2010; DIO, and 
references therein), the metallicity-HB index diagram indicates that 
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Figure 13. The metallicity-HB morphology diagram for 65 Galactic Glob- 
ular Clusters from DIO is presented. Solid circles and open squares repre- 
sent inner (Rg^ < 8 kpc) and outer {Rg^ > 8 kpc) halo GGCs, respectively. 
The solid line is the fiducial fit for the inner GGCs from DIO. And XI and 
And XIII are marked as blue open stars along with 10 other Local Group 
dSph galaxies (Harbeck et al. 2001; Sarajedini et al. 2002; Komiyama et al. 
2007). The Local Group dSph galaxies generally exhibit redder HBs than 
the GGCs at the same metal abundance. 



Table 3. HB parameters 



Object 


< V{HB) > 


{V-Dhb 


(V-Drgb 


A(V - /) 






median 


median 




And XI 


25.31±0.02 


0.559±0.070 


1.108±0.020 


0.549±0.074 


And XIII 


25.49±0.02 


0.572±0.109 


I.06I±0.013 


0.489±0.110 



the stellar populations in these Local Group dwarf galaxies are ~ 
1-3 Gyr younger than the GGCs in terms of their mean ages (DIO). 
However, unlike GGCs, the vast majority of which are single-aged 
systems, the dSph galaxies are usually multiple-generation systems 
with various stellar populations formed through extended star for- 
mation and complex chemical enrichment processes. Considering 
this fundamental difference between the Local Group dwarf galax- 
ies and GGCs, age itself is not likely to be enough to fully explain 
the diversity in the HB morphologies of these dwarf galaxies. The 
metallicity-HB morphology relation of the Local Group dSph satel- 
lites likely reflects a combination of etfects from multiple param- 
eters such as [a/Fe], helium content, mass loss, and environmen- 
tal conditions (e.g. photo-ionizing radiation from host galaxies) - 
in other words, their star formation histories. In fact, from Figure 
13, it would appear that dSphs with a lack of recent star formation 
events tend to show bluer HB morphologies (Sculptor, And XI, and 
And XIII), while the dSphs with extended and more recent star for- 
mation exhibit generally redder HBs. For example, according to 
Dolphin et al. (2005), over 90% of the stars in the Sculptor dSph 
formed more than 10 Gyr ago and subsequent star formation has 
largely ceased so that there are few signs of recent star formation 
events. As a result, the HB morphology of the Sculptor dSph re- 
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Figure 14. The correlation between A(V-I) and (B-R)/(B+V+R) index (see 
Figure 2 of DIG). The dashed line illustrates the least square fit of a cubic 
polynomial to the data. We used this relation to convert the (B-R)/(B+V+R) 
values of 8 other Local Group dSph galaxies obtained from Harbeck et al. 
(2001). 
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Figure 15. The measurement of (B-R)/(B+V+R) indices for And XI and 
And XIII. The blue and red HB stars are selected by the enclosed regions. 
The boundaries of the instability strip were obtained from the (V-I) colors 
of RRL stars. 



mains mostly blue and has not been significantly contaminated by 
intermediate aged RC stars compared to other dSphs. In addition, 
as mentioned in the previous section, both And XI and And XIII 
show no indications of intermediate or young stellar populations in 
their CMDs. As such, their CMDs are reminiscent of purely old, 
metal-poor GGCs. 

In order to probe this effect in more detail, we note that Dol- 
phin et al. (2005) find that for dSph and dwarf elliptical systems, 
the more intrinsically luminous galaxies generally exhibit extended 
star formation episodes, while the less luminous ones are more 
likely to have formed the bulk of their stellar populations over 10 
Gyr ago. This is graphically illustrated in the upper panel of Fig- 
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Figure 16. The upper panel shows the relation between the star formation 
history of dwarf galaxies over the past 10 Gyr (/iog) from the work of Dol- 
phin et al. (2005) and their integrated absolute magnitudes. For the lower 
panel, we measure the difference between each dwarf galaxy's HB index 
and that of the inner halo globular clusters in Figure 1 3 [A(Fit - Measured)] 
and plot that as a function of the galaxy's absolute magnitude. 



ure 16 where we have plotted the data in Table 1 of Dolphin et al. 
(2005). The /loo quantity represents the fraction of the stellar pop- 
ulation in each galaxy that has formed between the present and 10 
Gyr ago. As such, a low value of /iog means that there has been 
little recent star formation. The upper panel of Figure 16 suggests 
therefore that there is a correlation between absolute magnitude and 
the amount of recent star formation. 

To examine the validity of our assertion - that the degree of 
recent star formation is responsible for the location of the dwarf 
galaxies to the left of the inner halo GC locus in Figure 13, we 
plot the deviation of each galaxy from this locus as a function of 
the galaxy's My (which is a reflection of the degree of recent star 
formation) in the lower panel of Figure 16. The error bars represent 
the metallicity errors of each system propagated through to errors 
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in A(Fit - Measured) (cf. Dotter et al. 2010). The lower panel of 
Figure 16 reveals that the deviation of each dSph from the inner 
halo line in Figure 13 seems to depend on its My (i.e. the degree 
of recent star formation). Thus, we have affirmed our assertion that 
the deviation of the dSphs from the inner halo line is likely to be 
due to their more complex star formation histories as compared to 
the majority of Galactic GCs. 

The one possible outlier in the lower panel of Figure 16 is the 
Sculptor dSph, which seems to have had a fairly quiescent star for- 
mation history in contradistinction to that implied by its My- This 
could imply that star formation in Sculptor was quenched by an 
external process such as gas stripping via a dynamical interaction 
with the Milky Way or another more massive galaxy. 



8 DISCUSSION 

The luminosity-metallicity (L-M) relation of the Local Group 
dwarf satellites reflects the intrinsic properties of dwarf galaxies 
and contains significant implications for the mechanism of galaxy 
formation and the environmental conditions present during the 
early epochs of star formation (Harbeck et al 2005; Bovill & Ri- 
cotti 2009). 

The standard scenarios of galactic chemical enrichment (Lar- 
son 1974; Tinsley & Larson 1979; Dekel & Silk 1986) suggest that 
a L-M relation should be present for dwarf galaxies. The L-M re- 
lation shows that more luminous (i.e. more massive) galaxies are 
more metal-rich while less luminous (i.e. less massive) ones tend 
to be more metal-poor (Mateo 2008). This agrees very well with 
the theoretical prediction that more massive systems largely retain 
their interstellar media in the face of galactic winds triggered by su- 
pernova explosions so that more metals are trapped in these dwarf 
galaxies (Saviane et al. 2008). The L-M relation of the Local Group 
dwarf galaxies has also revealed the fundamental differences and 
possible correlations between dwarf irregular (dirr) and dSph sys- 
tems (Skillman & Bender 1995; Richer et al. 1998; Mateo 1998; 
G03). From their comparative study of 40 nearby dwarf galaxies, 
G03 noted that the old stellar populations of dirrs are systemati- 
cally more metal-poor as compared with those of dSphs at the same 
luminosity (i.e. present-day dIrrs are at least 10 times brighter in 
luminosity than present-day dSphs at a given metallicity). This in- 
dicates that the gas deficient dSphs in which recent star formation 
has been largely dormant, should be more effectively enriched than 
dIrrs, which currently exhibit on-going star formation. Based on the 
present-day luminosity, H I gas content, and the modest rates of on- 
going star formation of the dIrrs, G03 concluded that most of the 
present-day dIrrs in the Local Group, cannot fade enough to evolve 
into dSphs over a Hubble time. Hence, the dIrrs are not likely pre- 
cursors of the dSphs; furthermore, these two classes of galaxies are 
fundamentally different stellar systems that have experienced dif- 
ferent paths in their evolution. 

With our newly determined metallicities from the RRL peri- 
ods, we have plotted And XI and And XIII on the L-M relation 
of the Local Group dwarf galaxies in Figure 17. The absolute V 
magnitudes of And XI and And XIII were adopted from the work 
of CIO. From among the 40 nearby galaxies in the compilation of 
G03, we have included 34 Local Group dwarf satellites within 1.5 
Mpc of the Milky Way (MW). The L-M relation for 8 recently 
discovered ultra-faint dwarf (uFd) satellites of the MW were ob- 
tained from the recent work of Norris et al. (2010). The ordinate 
of this plot represents the mean metallicity of the old stellar pop- 
ulations in each galaxy while the abscissa is the integrated abso- 
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Figure 17. Luminosity-metallicity (L-M) relation for Local Group dwarf 
galaxies. The data were obtained from the compilation of Grebel, Gallagher, 
& Harbeck (2003) for the canonical dwarf satellites {My < -8). For the uFd 
systems, the data come from the recent work of Norris et af (2010). And XI 
and And Xni are shown as marked boxes. The metallicity-Mi/ relation ap- 
pears to be very well established from the brightest dE, M32, to the faintest 

MW satellite. Segue 1. There is a noticeable gap around My = -6 8. 

And XI and And XIII seem to fill this gap but they exhibit relatively higher 
metallicities than the mean L-M relation at a given absolute V magnitude. 
The typical error in the mean metallicity is about ~ 0.3 - 0.4 dex. 



lute V magnitude of each dwarf galaxy. To better understand this 
plot, a different symbol was used for each class of dwarf galaxy 
(solid circles (dIrr); open circles (dSph); open triangles (dE); open 
stars (dIrr/dSph, transition type); solid squares (ultra faint dwarf)). 
From Figure 17, we see that even given a typical metallicity error 
of cr^[Pc/H]> ~ ±0.4 dex, a well established L-M relationship exists 
from the brightest dIrr and dE to the faintest dwarf satellite. Segue 
1. However, there is a noticable gap between the faint end of the 
canonical dwarf galaxies and the bright end of the newly discovered 
ultra-faint dwarf satellites. We also observe the well known fea- 
tures in the L-M relationship that dSphs are more metal-rich than 
the dIrrs at the same absolute V magnitude, while the transition 
type dwarf galaxies (dIrr/dSph : Phoenix, DDO 210, LGS 3, and 
Antlia) show similar characteristics to the dSphs (Skillman & Ben- 
der 1995; Richer et al. 1998; Mateo 1998, and G03). Our estimates 
of the metallicity (< [Fe/H] >A„dxj=-l-^5; < [Fe/H] >A„dxin=- 
1.74) tend to locate And XI and And XIII slightly off from the 
average L-M relation as compared with the previous work of CIO 
(< [Fe/H >=-2.0 for both And XI and And XIII). However both 
galaxies still appear to follow the general trend of this L-M relation 
within the errors of the metallicity determinations. 

The observed properties of the newly discovered uFd galaxies, 
such as size, surface brightness, mass to light ratio, and L-M rela- 
tion can be explained within the context of "the tidal disruption sce- 
nario" of dSph formation which suggests that the present-day faint 
galaxies were once much more massive but lost significant mass 
via tidal interactions with their giant host galaxies (e.g. Ursa Minor 
dSph : Martinez-Delgado et al. 2001; Mayer et al. 2001). However, 
this tidal scenario cannot account for all of the properties of uFd 
galaxies, especially the very metal-poor nature of some of the uFd 
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Table 4. HB index and other physical parameters of sample dSph galaxies. 



Object 


(B-R)/(B+V+R) 


A{V - 7) 


[Fe/H] 


[M/H] 


Mv 


Sculptor 


0.06 


0.470 


-1.5 ±0.5 


-1.40 


-9.8 


Sextans 


-0.37 


0.270 


-1.9 ±0.4 


-1.80 


-9.5 


Tucana 


-0.20 


0.312 


-1.7 ±0.2 


-1.60 


-9.6 


And I 


-0.80 


0.198 


-1.4 ±0.2 


-1.30 


-11.8 


And II 


-0.70 


0.212 


-1.5 ±0.3 


-1.40 


-11.8 


And III 


-0.67 


0.217 


-1.7 ±0.2 


-1.60 


-10.2 


And V 


-0.62 


0.224 


-1.9 ±0.1 


-1.80 


-9.1 


And VI 


-0.70 


0.212 


-1.7 ±0.2 


-1.60 


-11.3 


And XI 


0.25 


0.549 


-1.75 ±0.12 


-1.65 


-7.3 


And XIII 


0.13 


0.489 


-1.74 ±0.12 


-1.64 


-6.9 



galaxies. If the dSphs formed in massive halos, the deep gravita- 
tional potential well of these halos should have effectively retained 
their metals. Hence, there should be a lower limit to the metallic- 
ity of the dSph galaxies. The continuously decreasing metallicity in 
the L-M relation from the bright dirrs and/or dEs to the uFd satel- 
lites may imply that the uFd galaxies were originally less massive 
stellar systems formed in low mass dark matter halos with fewer 
metals in which most of stars formed before the end of reionization 
(6 < z < 20; Bovill & Ricotti 2009; Walker et al. 2009; Salvadori 
& Ferrara 2009). In terms of the absolute V magnitude, And XI 
(Mv=-6.9) and And XIII (Mv=-6.7) appear to reside in the uFd 
regime (CIO). As mentioned above, their metallicities seem to be 
relatively higher than the average L-M relation at a given luminos- 
ity. This opens up the possibility that these two faint dSph may 
be tidal remnants from interactions with their giant host, M31 or 
prereionization fossils in which metallicity was enriched by a late 
phase gas accretion and star formation event well after reioniza- 
tion at a redshift z< 1-2 (Ricotti 2009). However, considering the 
likelihood that the absolute V magnitudes of And XI and And XIII 
are lower limits (M06), we cannot ignore the possibility that these 
faint dSph could be canonical dSph galaxies. Therefore, the verdict 
about the true origin of And XI and And XIII is still tentative. 



9 SUMMARY AND CONCLUSIONS 

We present a comprehensive study of the stellar populations in 
two faint M31 dwarf satellites. And XI and And XIII, using deep 
archival images from HST/WFPC2. Based on the analysis of the 
HB morphology and the properties of the RRL stars, we obtain the 
following results : 

1. Based on the appearance of their CMDs, And XI and And 
XIII appear to be purely old and metal-poor stellar systems and ap- 
pear to have experienced relatively simple star formation histories 
as compared with their more luminous counterparts in the Local 
Group. 

2. We have discovered and characterized RRL populations in 
both And XI (N„i, = 10, N, = 5) and And XIII (N„i, = 8,N, = 1) 
using our new template light curve fitting routine (RRFIT). The 
Period- Amplitude relations of the RRL stars in these two M31 
dwarfs appear to follow the P-A trend of the RRL population found 
in the Brown field (B04). The lack of AC populations reinforces 
our assertion that both And XI and And XIII are purely old stellar 
systems 

3. The metallicities of the RRab stars were calculated via the 
P-A-[Fe/H] relationship of Alcock et al. (2000). The metallici- 



ties obtained in this way {[Fe/H]A„jxi = -1.75; [Fe/H]A„jxiii = 
-1.74) are consistent with the values calculated from the RGB 
slopes indicating that our measurements are not significantly af- 
fected by the evolution of the RRLs away from the zero age hori- 
zontal branch. 

4. The distance to each galaxy was determined using 
the absolute V magnitudes of the RRab stars. We find (m - 
M)o.v/=24.33±0.05 for And XI and (m-M)o,v=24.62±0.05 for And 
XIII. Our results agree very well with the measurements from pre- 
vious studies (M06, CIO). 

5. The HB morphologies of And XI and And XIII were evalu- 
ated with a new metric called the A(V-I) index (DIO); this indicates 
that both of these faint galaxies are purely old systems. Our com- 
parative analysis of the HB morphology also revealed that the frac- 
tion of the ancient stellar populations (> 10 Gyr) in a Local Group 
dSph galaxy can be estimated by its HB morpholgy index. 

6. The locations of And XI and And XIII in the luminosity- 
metallicity plane follow the general trend for the Local Group 
dwarf galaxies and appear to fill the gap between the faint end of 
the canonical dSphs and the bright end of the newly discovered uFd 
systems. The metallicities of And XI and And XIII seem to be rel- 
atively higher than the average metallicity of uFd systems at the 
same luminosity. This may be an indication that these dwarfs were 
more massive in the past having experienced significant mass loss 
via tidal interactions with their giant host M31. Another possibility 
is that these systems are reionization fossils in which metallicity 
was enriched due to late phase gas accretion and star formation 
well after reionization (z ~ 1-2). A third possibility is simply that 
there are significant uncertainties in the integrated absolute V mag- 
nitudes of And XI and And XIII and that the true values are brighter 
than the measured ones. 
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APPENDIX A: THE SPECIAL CASE OF AND II 

In examining the appearance of Eig. 14, it seems that And II is 
the only M3 1 dwarf that exhibits a significant period spread among 
its ab-type RRLs in the P-A diagram. This was noted by P04 as 
well who suggested a correlation between the ab-type RRL period 
spread and the significant metallicity dispersion seen among the 
RGB stars in the And II CMD. We would like to investigate the 
bimodality of this period distribution and to correlate the inferred 
metallicity of the two RRL populations to see if it is consistent with 
the appearance and location of the And II RGB. 

Figure Al shows the P-A diagram for And II wherein we have 
fitted a straight line to the ab-type RRL relation. For each such star, 
we then calculate the offset between the measured period and the 
predicted period (ALog P) from this least squares fit, which is plot- 
ted as a solid line in Fig. Al. The inset to Figure Al shows the 
histogram of ALog P values generated in this way. Fitting a single 
Gaussian to this histogram (dashed line in Figure Al inset) pro- 
duces a reduced-;j'^ value of ~14 which suggests that there is a less 
than 1% chance that the distribution of ALog P is unimodal. In con- 
trast, fitting the ALog P histogram with a bimodal Gaussian yields 
a reduced-;('- of ~0.3, suggesting good agreement between the two 
distributions. The peaks of the two Gaussians differ by 0.100 ± 
0.006 dex in Log P, which corresponds to a metallicity difference 
of -0.90 ± 0.05 dex in [Fe/H] according to the Alcock et al. (2000) 
equation quoted above. 

Turning now to the appearance of the And II RGB, in Figure 
A2 we show the distance and reddening corrected CMD of And II 
along with the standard globular cluster RGBs from Sarajedini & 
Layden (1997). The And II photometry represents our reduction of 
the program frames taken for GO-6514 using the same procedure 
and techniques as described in Sec. 2 above. The sequences for 
M15, NGC 6752, NGC 1851, and 47 Tuc are plotted as grey lines. 
According to Zinn & West (1984), these have metallicities of -2.17, 
-1.54, -1.29, and -0.71, respectively. Our aim is to compare the 
metallicities from the RGB with the RRL values from the Alcock et 
al. (2000) relation, which is also based on the Zinn & West (1984) 
abundance scale. The GGCs have been corrected for distance using 
the V(HB) and E(B-V) values from Sarajedini & Layden (1997) 
along with the RRL luminosity-metallicity relation from Chaboyer 
(1993, see below). For And II, we adopt V(HB) = 24.93 and E(B- 
V) = 0.06 from values from Da Costa et al. (2000). 

The appearance of the And II CMD suggests the presents of 
a bifurcation in its RGB brightward of My 1.5 with a break es- 
sentially midway between the NGC 6752 and NGC 1851 RGBs. 
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Figure Al. The Period-Amplitude diagram for And 11. The solid hne rep- 
resents a linear least square fit to the ab-type RRL relation. The inset his- 
togram illustrates the distiibution of A log P values which is the offset of 
the measured period for each RRab star from the least square fit. A single 
Gaussian fit to this histogram is shown as a dashed line, while a bimodal 
Gaussian is described by a solid line. 

As such, we have interpolated between these two RGBs and con- 
structed a mean RGB represented by the sohd black line in Fig. A2. 
The solid line in the inset to Fig. A2 shows the histogram of color 
differences between this mean RGB sequence and each of the 134 
RGB stars brighter than Mv=-1.5. The dashed histogram in the in- 
set is constructed using the original photometric data consisting of 
74 RGB stars from Da Costa et al. (2000) showing good agreement 
with our result. A bimodal Gaussian fit to the solid line in the Fig. 
A2 inset suggests a difference of 5[A(B-V)] = 0.352 ±0.016 in the 
peaks of the two distributions. Using Equation (4) of Sarajedini & 
Layden (1997), this corresponds to a metallicity difference of ~0.70 
± 0.03 dex in [Fe/H]. In terms of the random errors, this value is 
3.4cr different than the metallicity difference inferred from the ab- 
type RRLs above, but given the fact that the systematic errors are 
likely to be larger than the random ones, the metallicity difference 
inferred from the RGB stars and the ab-type RRLs are consistent 
with each other. 

This paper has been typeset from a TjX/ KTjX file prepared by the 
author. 
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Figure A2. The distance and reddening corrected CMD of And II. The RGB 
sequences for Ml 5, NGC6752, NGC1851, and 47 Tuc obtained from Sara- 
jedini & Layden (1997) are illustrated as grey fines. The solid black line 
represents mid-points between the NGC6752 and NGC 1851 RGBs, which 
is coincident with the bifurcation in the RGB of And II brightward of My 

1.5. The solid histogram in the inset figure illustrates the distribution of 

(B-V) color dilferences between this mean RGB sequence (solid black line) 
and And II RGB stars brighter than Mi/=-1.5. The dashed histogram in the 
inset is constructed using the original And II data as published by Da Costa 
et al. (2000). The solid line in the inset represents a bimodal Gaussian fit to 
the solid histogram. 



